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Abstract. - The heavy-fermion metal CePdi-xRhx evolves from ferromagnetism at x = to a 
non- magnetic state at some critical concentration Xc- Utilizing the quasiparticle picture and the 
concept of fermion condensation quantum phase transition (FCQPT), we address the question 
about non-Fermi liquid (NFL) behavior of ferromagnet CePdi_xRhx and show that it coincides 
with that of both antiferromagnet YbRh2(Sio.95Geo.05)2 and paramagnet CeRu2Si2 and CeNi2Ge2. 
We conclude that the NFL behavior being independent of the peculiarities of specific alloy, is 
universal, while numerous quantum critical points assumed to be responsible for the NFL behavior 
of different HF metals can be well reduced to the only quantum critical point related to FCQPT. 
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' The nature of quantum criticality determining the non- 
f'ermi liquid (NFL) behavior observed in heavy-fermion 
(HF) metals is everyday topic of the physics of correlated 
Electrons. A quantum critical point (QCP) can arise by 
suppressing the transition temperature Tc of a ferromag- 
i).etic (FM) (or antifcrromagnctic (AFM)) phase to zero 
by tuning some control parameter ^ other than temper- 
ature, such as pressure P, magnetic field B, or doping 
X as it takes place in the case of the HF ferromagnet 
CePdi_xRhx [1,2] or the HF metal CeIn3_xSnx [3]. The 
NFL behavior around QCPs manifests itself in various 
anomalies. One of them is power in T variations of the 
specific heat C{T), thermal expansion a{T), magnetic sus- 
ceptibifity x{T) etc. 

It is widely believed that the NFL behavior is deter- 
mined by quantum phase transitions which occur at the 
corresponding QCP's. According to this concept, NFL be- 
havior in this case is due to the presence of thermal and 
quantum fluctuations suppressing quasiparticles [4-6] so 
that the quantum criticality in these systems can be de- 
scribed by conventional theory related to a spin-density- 
wave instability [7] or scenarios where the heavy electrons 
localize at magnetic QCP's, for example, due to a dcstruc- 
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tion of the Kondo resonance [8] . Unfortunately, up to now 
it was not possible to describe all available experimental 
facts related to the NFL behavior within a single theory 
based on the above scenarios. 

Measurements performed on the three dimensional FM 
CePdi_xRhx show that around some concentration x — 
Xc — 0.87 — 0.9 the suppression of the FM phase oc- 
curs, so that this alloy is tuned from ferromagnetism at 
X = to a non-magnetic state at QCP with the crit- 
ical concentration [1,2]. At a; = Xc, measurements 
on CePdi_xRhx show that the electronic contribution to 
the specific heat C(T) and the thermal expansion coeffi- 
cient a{T) behave as C(T) (x a(T) oc ^/T [1,9]. At the 
concentrations x < Xc, C(T)/T shows a peak at some 
temperature Tmax, while under the application of mag- 
netic field Tniax shifts to higher values [2]. Above dis- 
cussed scenarios for NFL behavior [6-8] imply that its 
details would in particular depend on system's magnetic 
ground state. Namely, within these scenarios, one can as- 
sume that the NFL peculiarities of CePdi_xRhx are to 
be different from those of either CeNi2Ge2 and CeRu2Si2 
exhibiting a paramagnetic ground state [10, 11] or from 
those of AFM cubic HF metal Celng.xSux [3] and HF 
metal YbRh2(Sio.95Geo.o5)2 exhibiting (in measurements 
of C{T)/T ) a weak AFM ordering at T < 20 mK [12]. 
On the other hand, the measurements of x{T^) have shown 
that the quantum critical fluctuations in this metal have a 
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strong FM component and thus are unique among all other 
quantum critical HF systems [13]. Obviously the critical 
fluctuations taking place at QCPs in the different HF met- 
als are different so that it may seem that we cannot have a 
universal behavior in these metals. Also, the above tradi- 
tional scenarios have no grounds to consider these QCPs 
as different manifestations of some single QCP. Moreover, 
the behavior of C[T)/T in YbRh2(Sio.95Geo.o5)2 is formed 
by AFM fluctuations while that of x{T) is determined 
by FM ones. The distinctive features of FM, AFM and 
paramagnetic systems suggest the intrinsic differences in 
their QCPs resulting in the diversity of their thermody- 
namic properties. Existing theories corroborate this point 
of view, they predict that magnetic and thermal proper- 
ties of CePdi_xRhx [1,2,4-6,14] should differ from those of 
YbRh2(Sio.95Geo.o5)2 since the latter substance is suppose 
to combine FM and AFM orders. 

Below we shall see that NFL properties of the function 
C{T)/T in CePdi_xRhx coincide with those of x{T) in 
CeRu2Si2 and YbRh2(Sio.95Geo.o5)2 as well as with those 
of C(r)/T in YbRh2(Sio.95Geo.o5)2- Also, the NFL be- 
havior of a{T) in CePdi_xRhx coincides with that of a{T) 
in HF metals CeNi2Ge2 and CeIn3_xSnx. The observed 
power laws and universal behavior of C{T) and a{T) in 
CePdi_xRhx can be hardly accounted for within the above 
scenarios when quasiparticles are suppressed, for there is 
no reason to expect that C(T), x(T), a{T) and other ther- 
modynamic quantities are affected by the fluctuations or 
localization in a correlated fashion. 

It might be possible to explain this universal behavior by 
Landau Fermi liquid (LFL) theory based on the existence 
of quasiparticles since C{T)/T cx a{T) cx x{T) cx M* 
where M* is the effective mass. Unfortunately, the effec- 
tive mass of conventional Landau quasiparticles is tem- 
perature, magnetic field, pressure etc. independent [15] 
and this fact contradicts to the measurements on HF met- 
als. On the other hand, when the electronic system of 
HF metals undergoes the fermion condensation quantum 
phase transition (FCQPT), the fluctuations are strongly 
suppressed and cannot destroy the quasiparticles which 
survive down to the lowest temperatures [16-19]. In con- 
trast to the conventional M* , the effective mass of these 
quasiparticles strongly depends on T, x, B etc. so that we 
have every reason to suggest that they are indeed respon- 
sible for the universal behavior observed in HF metals. 
We note that the direct observations of quasiparticles in 
CeCoIns have been reported recently [20]. 

In this Letter we show that the NFL properties of HF 
metals coincide regardless of their magnetic ground state 
properties. Namely, the NFL features observed in FM 
CePdi_xRhx, in cubic AFM CeIn3_xSnx, in paramagnets 
CeNi2Ge2 and CeRu2Si2 and in YbRh2(Sio. 95660. 05)2 dis- 
playing both AFM and FM fluctuations, coincide. Our 
main conclusion is that observed universal behavior is in- 
dependent of the peculiarities of the given alloy such as 
its lattice structure, magnetic ground state, dimensional- 
ity etc. so that numerous previously introduced QCPs can 



be substituted by the only QCP related to FCQPT. 

The schematic phase diagram of the HF metals under 
consideration is reported fig. 1. We show two LFL re- 
gions (left one being paramagnet (PM) or having long- 
range magnetic order and right one corresponds to reen- 
trant LFL phase induced by a magnetic field), separated 
by NFL one. The control parameter C, (see also above) 
can be pressure P, magnetic field B, or doping x. The 
variation of C drives the system from LFL region to NFL 
one and then again to LFL. The caption "Magnetic field 
induced LFL" means that only magnetic field can gener- 
ate the reentrant LFL phase. If C, is not a magnetic field, 
the right LFL-NFL boundary lies on the abscissa axis. 

To study the universal low temperature features of HF 
metals, we use the model of homogeneous heavy-electron 
liquid with the effective mass Af *(T, B, p), where the num- 
ber density p ~ pp/Sir'^, and is the Fermi momen- 
tum [15]. This permits to avoid complications associated 
with the crystalline anisotropy of solids [17]. To describe 
the effective mass M* (T, B) as a function of temperature 
and applied magnetic field B, when the heavy-electron sys- 
tem evolves from the LFL state, we use the Landau equa- 
tion relating the effective mass Af *(T, B) to the bare mass 
M and Landau interaction amplitude P(pi,p2,p) [15] 
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where n(p, T, R) is the quasiparticle distribution function 
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R = hbB/T. Here ^ = e(p, T) — fi, jis is the Bohr mag- 
neton, e(p,T) is the single-particle energy and p, stands 
for a chemical potential. 

We first consider the case when at T ^ the hcavy- 
clcctron liquid behaves as LFL and is located on the Fermi- 
liquid (FL) side of FCQPT (see Ref. [18] for details). Since 
e{p = pf) = p at B ^ 0, we see from eq. (3) that 
n{p,T,B) 6{pf —p), 6{p) is the step function. In this 
case eq. (1) reads [15,21] 
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Here iVg is the density of states of a free electron gas, 
F^{pf,Pf,p) is the p-wave component of Landau ampli- 
tude. LFL theory implies that the amplitude can be rep- 
resented as a function of p only, F^{pf,Pf, p) = F^{p)- 
We assume that at p ^ ppc, F^{p) achieves some value 
where the denominator tends to zero and find from eq. (4) 
that the effective mass diverges as [22,23] 
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where A, Ai are constants and ppc is QCP of FCQPT. 
Assuming that the control parameter ^ is represented by 
X and Xc corresponds to ppc we obtain {Cfc ~ C) /Cfc = 
{xc - x)/xc ~ (pfc - p)/pfc, while at C > Cfc the system 
is on the fermion condensation (FC) side of FCQPT [18]. 

Now we consider the temperature behavior of the ef- 
fective mass M*{T) in a zero magnetic field. Upon us- 
ing eq. (4) and introducing the function (5n(p,T) = 
n(p, T) — 9{pf — p), eq. (1) takes the form 
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We integrate the second term on the right hand side of eq. 
(6) over the angular variable 17, use the notation 



Fi {pF,p,p) = MpF j Pf 



aF(pF,p,p) 



apF (2^)3' 

and substitute the variable p hy z = ^{p) /T . Since in HF 
metals the band is flat and narrow, we use the approxi- 
mation ^{p) ~ pp{p — pf)/M* (T) and with respect to eq. 
(6) finally obtain 
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Here /3 = TM*{T)/p\ and f{z) = Fi{pF,z,p). The mo- 
mentum Pf is defined from the relation £{pf) = P- 

To investigate the low temperature behavior of M*{T), 
we evaluate the integral (8). Going beyond the usual ap- 
proximation [24] , we may obtain following final result 

/3/(0)ln{l + cxp(-l//3)} 



M*(T) 
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where Ai and A2 are constants of order unity. Here 
the logarithmic term is the result of an effective summa- 
tion of the main nonanalytic (at T 0) contributions, 
proportional to exp(— To analyze eq. (9), we first 
assume that /? <C 1. Then, omitting terms of the order 
of exp(-l//3), we obtain that at T < p\,/M*{p) 
the sum on the right hand side represents a T^-correction 
to M*{p) and the system demonstrates the LFL behav- 
ior [25]. At higher temperatures, the system enters a 
transition regime when the effective mass reaches its max- 
imal value M'^j at some temperature Tm- It can be easily 
checked that the terms proportional to 0^ and in eq. 
(9) are "responsible" for the maximum. The normalized 
effective mass M]^(T) = M*{T)/Mlj as a function of nor- 
malized temperature T/v = T/Ti\i is reported in the inset 
to fig. 1, showing several regimes. At T/v <C 1, the LFL 
regime with almost constant effective mass, occurs. At 
Tff ~ 1 it gives place to the transition region. At ele- 
vated temperatures when M/M*{p) <^ f3^, eq. (9) reads 
M/M*{T) (X T'^M*{Tf, giving [25,26] 

(10) 
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Fig. 1: Schematic phase diagram of the systems under consid- 
eration. Control parameter C, represents doping x, magnetic 
field B, pressure P etc. ("fc denotes point at which the ef- 
fective mass diverges. If C, is not a magnetic field, then the 
right boundary line NFL-LFL lies on the abscissa axis. Inset 
- normalized effective mass M^{T) = M*{T)/Mm {Ml, is its 
maximal value at T = Tm) versus the normalized temperature 
Tjv = T/Tm- Several regions are shown. First goes the LFL 
regime {Mlf{T) ~ const) at Tjv ^ 1, then transition regime 
(shaded area) where M^jiT) reaches its maximum. At elevated 
temperatures T^^/'^ regime given by eq. (10) occurs followed 
by T^^l"^ behavior, see eq. (11). 



Numerical calculations based on cqs. (8) and (9) show 
that at rising temperatures the linear term cx /3 gives the 
main contribution and leads to new regime when cq. (9) 
reads M/M*{T) cx /? yielding 



M*(T) cx r"^/2. 



(11) 



Note, that "rising temperatures" are still sufficiently low 
for the expansion of integrals in eq. (8) in powers of (3 to be 
valid. In the inset to fig. 1 both T^"^/^ and T^^l"^ regimes 
are marked as NFL ones since the effective mass depends 
strongly on temperature, which is not the case for the 
transition region. If the system is located at the FCQPT 
critical point, it follows from cq. (5) that Af*(pFc) ^ 00 
and Tp making the LFL region vanish, while the 
behavior of the effective mass at finite temperatures is 
given by cq. (11) [27]. The application of magnetic field 
restores the LFL behavior and at T = the effective mass 
depends on B as [26,28] 



M*{B) cx (B-B,o) 



-2/3 
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where Bco is the critical magnetic field driving both HF 
metal to its magnetic field tuned QCP and corresponding 
Neel temperature toward T = 0. In some cases Bco = 
0. For example, the HF metal CeRu2Si2 is characterized 
by Bco = and shows neither evidence of the magnetic 
ordering or superconductivity nor the LFL behavior down 
to the lowest temperatures [11]. In our simple model Bco is 



p-3 



V.R. Shaginyan et al. 



taken as a parameter. At elevated temperatures and fixed 
magnetic field, the effective mass depends on temperature 
as in the case when the system is placed on the FL side 
in accordance with eqs. (10) and (11) [25,29]. Since the 
magnetic field enters eq. (1) as the ratio R = fisB/T, at 
T/v < 1 the behavior of the effective mass can be described 
by a simple function 



M*{B,T) l + Cl7^2 



M*{B) l + C27^8/3' 



(13) 



which represents an approximation to solutions of cq. (1) 
that agrees with eqs. (10) and (12). Here 7^ = T/[{B - 
^00)^*3] I ci and C2 are fitting parameters. As we have seen 
the effective mass reaches its maximal value M^j at some 
TZ = TZm and we again define a normalized effective mass 
as M^{T,B) = M*{T,B)/M*j. Taking into account eq. 
(13) and introducing the variable y = TZ/TZm wc obtain 
the function 



M*{B) 1 + C12/2 
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which describes a universal behavior of the effective mass 
M^{y) when the system transits from LFL regime to that 
described by eq. (11). At p < ppc, ^I*{p) is finite, see 
eq. (5). In this case the eq. (14) is valid at T/v ^ 1 
if M*{T,B)/M*{p) < 1 because the term l/M*{p) on 
the right hand side of eq. (6) is small and can be safely 
omitted [27]. As a result, the behavior of M'^{y) has to 
coincide with that of the normalized effective mass M^{T) 
displayed in the inset to fig. 1. 

The effective mass M*{T,B) can be measured in ex- 
periments on HP metals. For example, M*{T,B) oc 
C{T)/T cx a{T)/T and M*{T,B) cx XAciT) where 
XAc{T) is ac magnetic susceptibility. If the correspond- 
ing measurements are carried out at fixed magnetic field 
B (or at fixed both the concentration x and B) then, as it 
follows from cq. (13), the effective mass reaches the max- 
imum at some temperature Tm- Upon normalizing both 
the effective mass by its peak value at each field B and the 
temperature by Tm, we observe that all the curves merge 
into single one, given by eq. (14) thus demonstrating a 
scaling behavior. 

As it is seen from fig. 2, the behavior of the normalized 
ac susceptibility XAciv) = Xac{TITm-,B)/xac{'^,B) ^ 
M^{Tn) obtained in measurements on the HP param- 
agnet CeRu2Si2 [11] agrees with both the approxima- 
tion given by eq. (14) and the normalized specific heat 
{C{Tn)/Tn)/C{1) = Mn{Tn) obtained in measurements 
on the HF FM CePdi_xRhx [2]. It is also seen from fig. 
2, that at temperatures T/v < 3 , the curve given by eq. 
(14) agrees perfectly with the measurements on CeRu2Si2 
whose electronic system is placed at FCQPT [29], that 
is in fig. 1 at C,yc- As to the normalized specific heat 
(shown by downright triangles in fig. 2) measured on 
CePdi-xRhx with x = 0.8 [2], the scaling holds up to 
relatively high temperatures. This is because its elec- 
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Fig. 2: Normalized magnetic susceptibility xn{Tn,B) — 
Xac{T/Tm,B)/xac{1,B) = M*t,{TN) for CeRuaSia in mag- 
netic fields 0.20 mT (squares), 0.39 mT (upright triangles) and 
0.94 mT (circles) against normalized temperature Tat = T /Tm 
[11]. The susceptibility reaches its maximum xac{Tm , B) at 
T = Tm- The normalized specific heat {C{Tn)/Tn)/C{1) of 
the HF ferromagnet CePdi_xRhx with x = 0.8 versus Tn is 
shown by downright triangles [2]. Here Tm is the temperature 
at the peak of C{T)/T. The solid curve traces the universal 
behavior of the normalized effective mass determined by eq. 
(14), it is also shown in figs. 3, 4, 5 and 6. Parameters ci and 
C2 are adjusted for xn{Tn, B) at B — 0.94 mT. 



tronic system is located on the FL side and the deflec- 
tion {xc — x)/xc ~ (p — Pfc) / Pfc at x = 0.8 from the 
critical concentration Xc — 0.9 is relatively big, elongating 
the T~^/^ region [27]. On the other hand, at diminish- 
ing temperatures the scaling is ceased at relatively high 
temperatures as soon as the LFL behavior related to the 
deflection from Xr sets in. 
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Fig. 3: The normalized effective mass M^{Tm , B) extracted 
from the measurements of the specific heat on CePdi-xRhx 
with X = 0.8 [2]. AT B > 1 T, MIj{Tn) coincides with that of 
CeRu2Si2 (solid curve, see the caption to fig. 2). 

Now we consider the behavior of M^(T), extracted from 
measurements of the specific heat on CePdi_xRhx under 
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the application of magnetic field [2] and shown in fig. 3. 
It is seen from fig. 3 that at S > IT the value M'^ de- 
scribes the normalized specific heat almost perfectly, co- 
incides with that of CcRu2Si2 and is in accord with the 
universal behavior of the normalized effective mass given 
by eq. (14). Thus, we conclude that the thermodynamic 
properties of CePdi-xRhx with x = 0.8 are determined by 
quasiparticles rather than by the critical magnetic fiuctu- 
ations. On the other hand, one could expect the growth 
of the critical fluctuations contribution as x —^ Xc so that 
the behavior of the normalized effective mass would devi- 
ate from that given by eq. (14). 




0.1 1 10 

Normalized temperature 

Fig. 4: Same as in Fig. 3 but a: = 0.85 [2]. At B > 1 T, 
Mj^(rjv) demonstrates the universal behavior (solid curve, see 
the caption to fig. 2). 

In fig. 4, the effective mass AI'^{Tn) at fixed i?'s is 
shown. Since the curve shown by circles and extracted 
from measurements at iJ = does not exhibit any maxi- 
mum down to 0.08 K [2], we conclude that in this case x 
is very close to Xc and function M'^{Tpf) is approximately 
described by eq. (11), while the maximum is shifted to 
very low temperatures or even absent. As seen from fig. 
4, the application of magnetic field restores the universal 
behavior given by eq. (14). Again, this permits us to con- 
clude that thermodynamic properties of CePdi_xRhx with 
X = 0.85 are determined by quasiparticles rather than by 
the critical magnetic fluctuations. 

The thermal expansion coefficient a(T) is given by [24] 
a{T) ~ M*T/{p%K{p)). The compressibility K{p) is not 
expected to be singular at FCQPT and is approximately 
constant [30]. Taking into account eq. (11), we find that 
a(T) (X Vt and the specific heat C(T) = TM* cx VT. 
Measurements of the specific heat C{T) on CePdi_xRhx 
with X ~ 0.9 show a power-law temperature dependence. 
It is described by the expression C{T)/T = AT^'^ with 
q ~ 0.5 and A=const [1]. Hence, we conclude that the 
behavior of the effective mass given by eq. (11) agrees 
with experimental facts. Measurements of a{T)/T on 
both CePdi_xRhx with x = 0.9 [1] and CeNi2Ge2 [10] 
are shown in fig. 5. It is seen that the approxima- 
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Fig. 5: The normalized thermal expansion coefficient 
(Q(r]v)/rjv)/a(l) = M*i^{Tn) for CeNi2Ge2 [10] and for 
CePdi_xRhx with x = 0.90 [2] versus Tjv = T/Tm- Data 
obtained in measurements on CePdi_xRhx at B = are mul- 
tiplied by some factor to adjust them in one point to the data 
for CeNi2Ge2. Dashed line is a fit for the data shown by the 
circles and pentagons at B = and represented by the func- 
tion a(T) — C3\/T with cs being a fitting parameter. The solid 
curve traces the universal behavior of the normalized effective 
mass determined by eq. (14), see the caption to fig. 2. 

tion a{T) = c^\/T is in good agreement with the results 
of measurements of a{T) in CePdi_xRhx and CeNi2Ge2 
over two decades in T^r. We note that measurements on 
CeIn3_xSnx with x = 0.65 [3] demonstrate the same be- 
havior a{T) (X y/T (not shown in fig. 5). As a result, 
we suggest that CeIn3_xSnx with x = 0.65, CePdi_xRhx 
with X ~ 0.9, and CeNi2Ge2 are located at FCQPT (in 
fig. 1 at Cfg) and recollect that CePdi_xRhx is a three 
dimensional FM [1,2], CeNi2Ge2 exhibits a paramagnetic 
ground state [10] and CeIn3_xSnx is AFM cubic metal [3]. 

Mj^(T/v) extracted from measurements on the HF met- 
als YbRh2(Sio.95Geo.o5)2, CeRu2Si2, CePdi_xRhx and 
CeNi2Ge2 is reported in fig. 6. It is seen that the uni- 
versal behavior of the effective mass given by eq. (14) is 
in accord with experimental facts. YbRh2 (Sip. 95009.05)2 
is located on the FC side where the system demonstrates 
the NFL behavior down to lowest temperatures [27]. In 
that case, C (see fig. 1) is represented by B and Cfc = -Bco- 
In the LFL regime induced by the magnetic field, the ef- 
fective mass M*{B) oc {B - Bcoy^/^ and does not foh 
low eq. (12) [27,29]. As a result, the range of the scal- 
ing behavior in temperature shrinks to the transition and 
2^-2/3 regions, see inset to fig. 1. It is seen from fig. 6 
that A/^(T7v) shown by downright triangles and collected 
on the AFM phase of YbRh2(Sio.95Geo.o5)2 [12] coincides 
with that collected on the FM phase (shown by upright 
triangles) of YbRh2(Sio.95Geo.o5)2 [13]. We note that in 
the case of LFL theory the corresponding normalized ef- 
fective mass M^j^ ~ 1 is independent of both T and B. 

The peak temperatures Tmaxi where the maxima of 
C{T)/T, XAc{T) and a{T)/T occur, shift to higher vah 
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Fig. 6: The universal behavior of AI^{Tn), extracted 
from xac{T,B)/xac{Tm,B) for both YbRh2(Sio.95Geo.o5)2 
and CeRu2Si2 [11, 13], {C{T)/T)/{C{Tm)/Tm) for both 
YbRh2(Sio.95Geo.o5)2 and CePdi-xRhx with x = 0.80 [2,12], 
and (a{T)/T)/{a{TM)/TM) for CeNi2Ge2 [10]. All the mea- 
surements are performed under the application of magnetic 
field as shown in the insets. The solid curve gives the universal 
behavior of M'^ determined by eq.(14), see the caption to fig. 2. 




B(T) 

Fig. 7: The peak temperatures Tmax(-B), extracted from mea- 
surements of C/T and xac on YbRh2(Sio.95Geo.o5)2 [12, 13] 
and approximated by straight lines. The lines intersect at 
B ~ 0.03 T. 

ues with increase of the applied magnetic field. It fol- 
lows from eq. (14) that Tm oc {B — Bco)^b- In fig. 
7, 71nax(5) are shown for C/T and xac, measured on 
YbRh2(Sio.95Geo.o5)2- R is seen that both functions can 
be represented by straight lines intersecting at ~ 0.03 T. 
This observation [12, 13] as well as the measurements on 
CePdi-xRhx, CeNi2Ge2 and CeRu2Si2 demonstrate the 
same behavior [2,10,11]. 

In summary, we have shown, that bringing the differ- 
ent experimental data (like C{T)/T, Xac{T), a{T)/T etc) 
collected on different HF metals (YbRh2(Sio.95Geo.o5)2, 
CeRu2Si2, CePdi_xRhx, CeIn3_xSnx and CeNi2Ge2) to 
the above normalized form immediately reveals their uni- 
versal scaling behavior. This is because all above experi- 



mental quantities are indeed proportional to the normal- 
ized effective mass. Since the effective mass determines the 
thermodynamic properties, wc conclude that above alloys 
demonstrate the universal NFL thermodynamic behavior, 
independent of the details of the HF metals such as their 
lattice structure, magnetic ground state, dimensionality 
etc. This conclusion implies also that numerous QCPs as- 
sumed earlier to be responsible for the NFL behavior of 
different HF metals can be well reduced to a single QCP 
related to FCQPT. 

This work was supported in part by RFBR, project No. 
05-02-16085. 
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